Quantum entanglement: The unitary 8- vertex braid matrix with imaginary rapidity 
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We study quantum entanglements induced on product states by the action of 8-vertex braid 
matrices, rendered unitary with purely imaginary spectral parameters (rapidity). The unitarity is 
displayed via the "canonical factorization" of the coefficients of the projectors spanning the basis. 
This adds one more new facet to the famous and fascinating features of the 8-vertex model. The 
double periodicity and the analytic properties of the elliptic functions involved lead to a rich structure 
of the 3-tangle quantifying the entanglement. We thus explore the complex relationship between 
topological and quantum entanglement. 

PACS numbers: 03.65.Ud Entanglement and quantum nonlocality; 03.65.Ca Formalism; 03. 67. -a Quantum 
information; 03.67.Mn Entanglement measures, witnesses, and other characterizations 
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"Are topological and quantum entanglements re- 
lated?" This intriguing question is recently being studied 
from different angles. One approach was intiated by Ar- 
avind (H- Kauffman and Lomonaco [H pointed out that 
braid matrices (representing the third Reidmeister move 
(sf, fundamental in the topological study of knots and 
links) correspond to universal quantum gates, when they 
are also unitary. In previous studies unitary braid matri- 
ces were constructed explicitly for all dimensions [4] and 
applied to the study of quantum entanglements [5] . Here 
our starting point is the 8-vertex model |6j with braid 
matrix related to the Yang-Baxter one through a suitable 
permutation of elements, rendered unitary by a passage 
to imaginary rapidity {0 iO). The consequent uni- 
tarity is displayed transparently through "canonical fac- 
torization" [7| of the coefficients of the projectors. One 
now no longer has a statistical model with real, posi- 
tive Boltzmann weights but unitarity thus implemented, 
opens a new road (as will be shown below) to quan- 
tum entanglements. We first formulate such unitariza- 
tion {0 iO) in a general fashion and illustrate with 
the relatively simple 6- vertex case. Then we concentrate 
on the far more complex 8-vertex case, and study the 3- 
tangle Pj parametrized by sums of products of ratios of 
the g'-Pochhammer functions. 

The R {0) being a N'^ x N"^ matrix, acts on the base 
space Vn Vn spanned by the tensor product of TV- 
dimensional vectors Vn- Defining R12 {0) = R{6) 
In, R23 {0) = In^R{0) , where In is the N x N identity 
matrix, the corresponding braid operator is 



B = Ru {e)R23 {0 + e')Ri2 {e') 

= R23 {0')Ri2 {0 + 0')R2s (0). 



(1) 



The above Braid equation corresponds to the equivalence 
of knots related through the third Reidemeister move [3] . 



Ref. provides an useful introduction to the equivalent 
Yang-Baxter formalism. Of course, B acts on the base 
space Vat Vat Vat . Additionally, if the braid matrix R 
is also unitary, then it induces unitary transformations 
in Vn ^ Vn, and B in Vn ^ Vn Vn- It is crucial to note 
the essential point that a non-trivial unitary R induces 
non-local unitarj^transformations. Had it beenjhe case 
that R = Ri ^ R2, where Ri is acting on Vi , R2 on V2 , 
and R'm Fi V2 , then such an R would have been trivial 
from the point of braiding. Thus a non-trivial B induces 
a non-local transformations in Vn Vn Vn- 

The non-local unitary actions set the stage for quan- 
tum entanglements. It was shown [5] that B, acting on 
un- entangled product states of the general form 



N 



N 



N 



in Vn (8) Vn ^ Vn, can generate entanglements for certain 
choices. We had also studied entanglements generated 
by two different classes (real and complex) of B. The 
"3-tangles" and "2-tangles" characterizing such entan- 
glements were obtained explicitly in parametrized forms 
in terms of the parameters of B, and the variations with 
{0,9') were analysed. 

In another paper [7] , we had introduced the "canonical 
factorization" for R{0), which turns out to be very sig- 
nificant. Here we will exploit all these results and show 
that the simple passage {0 iO) is sufficient to provide 
unitarity under the following constraints: 



(i) R{0) = 



fii-O) 
Pi = In2 {= In <S> In) ■ 

(ii) (RiO)) 

trans 

= Ri0). 



Pi, where PiPj — SijPi, and 



2 



Thus, initially R (0) is real and symmetric, with a com- 
plete set of orthonormal projectors Pi as a basis. The 
domain of i depends on the class considered. The fac- 
torized form fi (0)/ fi {—0) of the coefficients in the first 
constraint might seem strongly restrictive, but in fact 
it was shown that it holds true for all well-known stan- 
dard cases, and new such cases were constructed [7], with 
the new term "canonical factorization" being introduced. 
One can easily check that a direct consequence of the 
constraints is R{6) R{—6) =1^1. After the passage 
{0 -> iO), since Pj's are real, one can easily show that 
{R{iO)y{R{ie)) =1^1, i.e., R{ie) is unitary. 

First, we demonstrate this formalism with the simpler 
case of the 6- vertex models. Following Ref. [7|], which 
contains an extensive classification of "canonical factor- 
ization" for all dimensions, we define the projectors: 
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and obtain for the ferroelectric case after {0 iO), with 
the real parameter 7, 



i?(z^)=Pl(+)+Pl(_) 



cosh 1(7 -i(9) ^ 
" cosh 1(7 + z6>) 
sinh |(7 — iO) 
sinh I (7 + zl9)^ 



2(-), 



(3) 



which evidently satisfies the unitarity constraint. 

Now, we proceed to the more complicated case of the 
8-vertex model. We again define the projectors as in 
Eq. (|2]). The coefficients are expressed [10] in terms of in- 
finite products ((7-Pochhammer functions), starting with 



(^;«)oo = ^(l-^«")■ 



(4) 



n>0 



Setting z = exp(^), the initial 8-vertex matrix is: 
R (0) = {a + d) Pi(+) + ia-d) Pi(_) + ic + b) P2(+) 

+ (C-6)P2(-), 

(5) 

where with supplementary real parameters p, q one ob- 
tains 0: 

{a±d) = ) ^^Y^-^ ^^ (6) 



We note that defining the numerators of the two equa- 
tions Eq. (j6j) and Eq. ([7j) as fi(±){z) and /2(±)(^) respec- 
tively, and using the fact that z = exp(^), we can express 



them as (a ± d) 



and (c ± b) 



which implies that the essential property of the coeffi- 
cients, "canonical factorization" , is preserved. 

After {0 iO) passage, we thus have (a ± (i) = 

-P^Kw,^ and (c±b) = I^^^lKl • Since the other pa- 
rameters are real, we can interpret the coefficients as new 
phases {a±d) = e*^(±) and {c±b) = e*^(±) , where the 
phase factors (^(±),^(±)) are complicated functions of 
ip^q^O). Note also that the coefficients under complex 

conjugation become (aid)* = )\'^}_\(f.i6^ — (a±<i)~ 

/2(±)(e--^) 



and (c ± b) 



= {c±b)~^ . Since the projec- 



/2(±)(e- 

tors are real and symmetric, we again have the unitarity 
{R {iO)y R (iO)) = I^I. This opens the door of a new do- 
main as a generator of quantum entanglements, as shown 
hereafter. 

Consider the base space that is 8-dimensional and 
spanned by the states |ei) |e2) (8) lea) = 1616263) , where 
6i = ±; z = 1, 2, 3. We will adopt a notation (|+) , |— )) 
(|1) , 1 1)) that generalizes smoothly to higher spins. The 
braid operator is 

B = B^ = {^ie)^h){h^R{i0^iO')){R{ie')^l2), (8) 
and the matrix 



R{ie) = 



a 

c b 

b c 

d 



(9) 



where {a±d) = e^^(±)(^), {c±b) = e^^(±)(^), and 
(^piq-^e'^;p) (Tphe'^^^p) 

V / 00 V / c 



[TP^Q (^piqe^^;p] 



i<^>(±)(0)^ g^e ±q 

q2e 2 ± q ^2 
{Tpq~^e'^;p)^ {^pqe-'^;p)^ 
{TPq~^e-'^;p)^ {^pqe'^;p)^ 



(10) 



One crucial fact is that R has non-zero elements only 
on the diagonal and the anti-diagonal. This effectively 
splits the base space into two 4-dimensional subspaces 
closed under the action of B. They are spanned respec- 
tively by (|111),|111),|111),|111)) and V^^) = 
(|111) , |111) |111) , |111)), corresponding to even and odd 
numbers of indices with bar. Moreover, for say 



(c±6) 



(qiz i±q izi) {Tpq-^z;p) ^ {Tpqz-^;p) 

[q-^z-^±q 22; 2 



B |111) = ai |111) |111) +71 |111) +(^1 |111) , (11) 

(tw"^^"^p)oo (tw^;p)oo * 

(7) B |111) = ai llll) + /3i |Ill) + 71 |lll) + (5i |111) , (12) 
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with the same coefficients (o^i, 71, (^i). More gener- 
ally, the symmetry of (j9j) ensures for 

B \ijk) = ci \ijk) + 02 \ijk) + C3 \ijk) + C4 \ijk) , (13) 

with z, j, /c = (1 or 1), the direct consequence 

B \ijk) = ci \ijk) + 02 \ijk) + C3 + 04 \ijk) . (14) 

The coefficients are conserved as above for (i^j^k) 
{i^j^k). Thus it is sufficient to evaluate the action of B 
on the subspace V(e) or V(^o) • 

To study the behavior of density matrcies and 3- 
tangles, we explicitly consider the action of B on the 
product state |1) |1) (g) |1) = |111), in the subspace 
V(e) 7 given by (pT]) . Some straightforward algebra gives 

Oil = /(+)/(+)/'(+) + f{-)f\-)9'\+) 

li=9{-)f\+)f'\-)^9{+)f\-)9'\-) 

h = /(-)/(+)/'(+) + f{+)f\-)9\+)^ (15) 

where we have used the phase factors ^ {±) and ^{±) to 
define 



/(±) 
/(±) 
/(±) 
9{±) 



The matrix {P12P12) has the following eigenstates 





-)(^) 
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gi<3&(+)(^) ±e^*^(- 
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-)(^') 



(16) 



such that {f-,f'-,f"){±) correspond respectively to ar- 
guments {0^0' ^{0 + 0')) with analogous notations for 
{g, g' ^ g")(^±). Starting with (pT]) and tracing out the third 
index, one obtains the density matrix 



P12 



aia\ ai5l 

PtYi 7171* 
al5i 616I 



(17) 



Defining 
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one then obtains the matrix 





P12P12 2 



<^i<^i^^i^i 









/^i/^i*7i7i* /5f/5i*7i* 
7?/^i*7r /^i/^r7i7i* 




(18) 
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(20) 



with the eigenvalues 4aialSlSl,0,4f3if3^j^j^,0, respec- 
tively. Implementing the results of [8] (as in [5]), the 3- 
t angle, invariant under permutations of the subsystems 
(1, 2, 3), is obtained as 



ri23 = l6{aialh(3ljljlSiSiy 



(21) 



Due to the unitarity of B (after 6 iO) in (pTj) 
a^al + Pipl + 7*7* + 5i5l = 1 and < ri23 < 1. As 
the parameters (p, ^, 0') vary, the 3-tangle ri23 varies 
in the domain [0,1]. The doubly periodic elliptic func- 
tions involved, expressed in terms of the g'-Pochhammer 
functions as in the ratios (fTQ|) , demand painstaking com- 
putations involving rather involved algebra. This is in- 
deed the real attraction of the unitarized 8-vertex case. 
^ One can study entirely analogously 
5(|111) , |111) , |111)) in the subspace V^e) im- 
plementing respectively the sets of coefficients 
{ai,Pi,ji,Si),i = 2,3,4 as given by 
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= f{+)9\+)f'\-) - 
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= 9{-)9 (+)/ (+) - 
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= 9{+)9\+)9'\+) - 
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£l II 

= 9(+)J (+)9 (-) - 
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(22) 

Figure 1 shows the rich structure with subtle varia- 
tions for ri23, by the action of B on the product states 
|111) , |111) , |111) , |111) in the subspace V^^)- We note 
that for (^ + ^') = 0, we have ri23 = 0, so that in the do- 
main (— 7r,7r) for both {0^0')^ there are diagonal lines of 
symmetry, with a line of zero value passing through the 
origin. One further notes that (^,^)± (^,<l>)± and 
hence (ri23 ^ ^23), for {p l/p,q l/q,0 -0). 
There are more intricate and subtle lines of symmetry 
as evident in figure 2, where we show the oscillations of 
ri23 between zero and unity, as a function of (p, q) for 
5 I 111). The results for the other subspace V(o), namely 
B{\lll) , |111) |111) , |111)), follows from the symmetry of 
V(e) and V(o) under the action of B as stated in (p!3|) and 



FIG. 1. Variations of the 3-t angle T123 as a function of {0^0')^ by the action of B on the product states |111) , , , 
in the subspace V(e)5 given by ()lip . The parameters p — 0.1, q — 0.5. 




FIG. 2. Variations (cross-sectional and top views) of the 3- 
tangle T123 as a function of (p, ^) for {0 = n/^^O' = 7r/6), by 
the action of B on the product state |111). 

(p!4|) . Combining these results one can then study the 
action of B on the general product state, namely 

B {(xi \l)+xi |1)) ® (yi \l) + yi |1)) ® {z^ |1) + |1))} , 

with some more straightforward algebra. 

A central feature of quantum entanglements induced 
by unitary braid operators are parametrizations of the 
quantifiers of entanglements. The present case is a rich 
and subtle example. Various entangled cases chosen with 
simple constant coefficients to assure certain interesting 
properties (such as \GHZ) and \W)) are thus seen to 
be imbedded in a continuum when approached via braid 
operators. Such a continuum provides a link between 
topological and quantum entanglements. Some entan- 
glements are inequivalent under locally unitary transfor- 



mations. Non-local unitary transformations, an intrinsic 
feature of braid matrices, provide precise explicit unifi- 
cations. We intend to study all such aspects more thor- 
oughly elsewhere. Another quite different perspective 
will be also shown to be provided by an entirely different 
class of braid matrices (6'(9(n)-type [ll^) again "unita- 
rized" by as above {6 iO). There will emerge a spin 
chain linked with a class of Temperley-Lieb algebra and 
display another possibility of our basic approach. More- 
over, not being restricted to the "diagonal-antidiagonal" 
form (illustrated above by the 8- vertex model) these uni- 
tarized SO{n) matrices will generate a broader class of 
entanglements. 
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